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p^. Abstract 

D ■ 

,-J^ I Motivated by recent observations of non-zero ^13 from the Daya Bay and RENO experiments, we 

^f\ • propose a renormalizable neutrino model with A4 discrete symmetry accounting for deviations from 

\/^ . the tri-bimaximal mixing pattern of neutrino mixing matrix indicated by neutrino oscillation data. 

00 ; 

T— I I In the model, the light neutrino masses can be generated by radiative corrections, and we show 

^. ■ 

^y^ ' how the light neutrino mass matrix can be diagonalized by the Pontecorvo-Maki-Nakagawa-Sakata 

O ■ 

p^ ■ mixing matrix whose entries are determined by the current neutrino data including the Daya Bay 

L^ , result. We show that the origin of the deviations from the TBM mixing is non-degeneracy of the 

^ I neutrino Yukawa coupling constants, and unremovable CP phases in the neutrino Yukawa matrix 

H ■ 
. . .' gi'^6 rise to both low energy CP violation measurable from neutrino oscillation and high energy CP 

violation. 
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I. INTRODUCTION 

Very recently, the Daya Bay Collaborations [l| announced 5.2cr observation of the non- 
zero mixing angle 6'i3 with the result given by sin^ 26*13 = 0.092 ± 0.016(stat) ± 0.005{syst) 
^. This result is in good agreement with the previous data from the T2K, MINOS and 
Double Chooz Collaborations 3], and the Daya Bay and RENO progresses have led us 
to accomplish the measurements of three mixing angles, 612,623 and 6*13 from three kinds 
of neutrino oscillation experiments. A combined analysis of the data coming from T2K, 
MINOS, Double Chooz and Daya Bay experiments shows ^ that 

sin^ 2^13 = 0.089 ± 0.016(0.047) , (1) 

or equivalently 

a _ O fifi°+0-77° (+2.14°) fr,\ 

at 1(7 (3cr) levels and that the hypothesis ^13 = is now rejected at a significance level 
higher than 6a. In addition to the measurement of the mixing angle ^13, the global fit of 
the neutrino mixing angles and mass-squared differences at la {3a) levels are given by 5|] 

f fi c;°+l-6° (+4.2°) ^TT 
o -OAnO+l-0° (+2.9°) n _ ^fi i °+3.5° (+7.0°) ^ _ J "^'^ -1.4° (-4.7°) ' ^^^ 

^12 - d4.U _09o („2 7.) , ^23 - 4b.l „4oo („75.) , ^^ - \ +1.70 (+4^0) 

[^ '"^ -1.5° (-5.5°) ' ^^ 
[^ ^.4U_QQ9 (_o 27) ' -^^ 

in which NH and IH stand for normal hierarchical neutrino spectrum and inverted one, 
respectively. The data in Eqs. fl2l[3l) strongly support that the tri-bimaximal (TBM) mixing 
pattern of the lepton mixing matrix 61] should be modified. There have been theoretical 
attempts to explain what cause the three mixing angles to be deviated from their TBM 

n 

values [7|]. 

Motivated by the measurements of ^13 from the Daya Bay and RENO experiments, we 
propose in this paper a renormalizable model with A4 discrete symmetry which gives rise 
to deviations from the TBM mixing indicated by the current neutrino data. In addition 



^ The RENO Collaboration also announced observation of the non-zero mixing angle 6*13 [2[ in consistent 
with the result from the Daya Bay Collaboration. 



to the leptons and the Higgs scalar of the standard model (SM), the model we porpose 
contains three right handed heavy Majorana neutrinos and several scalar fields which are 
electroweak singlets required to construct desirable forms of the letponic mass matrices. 
Although we introduce electroweak singlet heavy Majorana neutrinos, the usual seesaw 
mechanism does not operate because the scalar field involved in neutrino Yukawa terms can 
not get vacuum expectation value (VEV). However, as will be shown later, the light neutrino 
masses can be generated through loop corrections which is a kind of the so-called radiative 
seesaw mechanism 8|. In the paper, we will show how the light neutrino mass matrix 
generated through loop corrections can be diagonalized by the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) mixing matrix whose entries are determined by the current neutrino data. 
The origin of the deviations from TBM mixing in our model is non-degeneracy of the neutrino 
Yukawa coupling constants among three generations, which is different from other attempts 
to explain the deviations from the TBM mixing [7|. 

Since non-trivial Dirac CP phase can exist only when the mixing angle ^13 has non-zero 
value in the standard parametrization of the leptonic mixing matrix, the observations of 
non-zero 613 from the Daya Bay and RENO experiments shed light on the search for CP 
violation in the leptonic sector. We will show that unremovable CP phases in the neutrino 
Yukawa matrix are the origin of the low energy CP violation measurable from neutrino 
oscillation as well as high energy CP violation. Therefore, we can anticipate that there may 
exist some correlation between low energy CP violation and high energy CP violation. 

II. A MODEL WITH A4 SYMMETRY 

The model we consider is the standard model (SM), extended to contain three right- 
handed 5'f/(2)L-singlet Majorana neutrinos, N^. In addition to the usual SM Higgs doublet 
$, we newly introduce two scalar fields, x ^^^ Vi that are singlet and doublet under SU(2)l, 
respectively: 

<^={v^,v'f , X, V={v\v'f ■ (4) 

In order to account for the present neutrino oscillation data, we impose ^4 fiavor symme- 
try for leptons and scalars. In addition to A4 symmetry, we introduce extra auxiliary Z2 
symmetry so that a radiative seesaw at around TeV scale should operate. Here we recall 



that A4 is the symmetry group of the tetrahedron and the finite groups of the even permu- 



tation of four objects 



3- The group A4 has two generators 5* and T, satisfying the relation 



52 
by 



T^ = (ST)^ = 1. In the three-dimensional unitary representation, S and T are given 



/ 



S 



\ 



\ 



( 



T 



1 

1 

1 



\ 



(5) 



1 
0-10 

-ly ^ ^ 

The group A4 has four irreducible representations, one triplet 3 and three singlets 1,1', 1" 
with the multiplication rules 3(g)3 = 3s©3a©l©l'© 1", 1' © 1" = 1, 1' © 1' = 1" and 
1" © 1" = V . Let's denote two A^^ triplets as (ai, a2, 03) and (61, 62, ^s), then we have 



(a © 6)3, 
(a © 6)3, 
(a ©6)1 
(a ©6)1' 
[a © h)ri 



(02^3 + «3^2, 0361 + 0163, 0162 + 02^1) 
(0263 - ^3^2, %&! - ai&3, ai^2 - 02^1^ 
aibi + 02^2 + a-zh , 
ai6i + a;a2fe2 + ^^03^3 , 



ai^i + ix'^a2&2 + 1^036; 



'3 , 



(6) 



where u 



3J27r/3 



is a complex cubic-root of unity. The representations of the field content 



of the model under SU{2) x f/(l) x A4 x Z2 are summarized in Table-H]: 



TABLE L Representations of the fields under A4 x Z2 and SU{2)l x U{1)y- 



Field 



^ei -^ui ^T 



7 // 7// 



''_R'''R 



iV^ 



X 



$ 



^4 


1, 1', 1" 


1, 1', 1' 


Z2 


-^ 


-F 


L X C/(l)y 


(2,-1) 


(1,-2) 



(1,0) 



3 


1 


+ 


+ 


L,0) 


(2,1 



(2,1) 



With the field content and the symmetries specified in Table [H the relevant renormalizable 
Lagrangian for the neutrino and charged lepton sectors invariant under SU{2) x [/(I) x ^4 x Z2 



is given by 

- £Yuk = yiLe{fjNR)i + y2Lf,{fjNR)v + ylLr{fjNR)r' 

+ Vel,^ Ir + y^L^^ 4 + y,L,$ l'^ + h.c , (7) 

where fi = zr2?7* with the Pauh matrix T2. Here, Le,^,T and /]j' denote left handed lepton 
SU{2)l doublets and right handed lepton SU{2)l singlets, respectively. The higher dimen- 
sional operators {d > 5) driven by x ^.nd 77 fields are suppressed by a cutoff scale A which 
is a very high energy scale. Thus, their contributions are expected to be very small and 
we do not include them in this work. In the above Lagrangian, mass terms of the charged 
leptons are given by the diagonal form because the Higgs scalar $ and the charged lepton 
fields are assigned to be A4 singlet. The heavy neutrinos Njn acquire a bare mass M as well 
as a mass induced by a vacuum of electroweak singlet scalar x assigned to be A^ triplet. 
While the standard Higgs scalar <l>° gets a VEV v = (2a/2G'f)"^/^ = 174 GeV, the neutral 
component of scalar doublet rj would not acquire a nontrivial VEV because 77 has odd parity 
of Z2 as assigned in Table |T] and the auxiliary Z2 symmetry is exactly conserved even after 
electroweak symmetry breaking ; 

(r/°) = 0, (z = l,2,3), ($0) = ,;^0. (8) 

Therefore, the neutral component of scalar doublet 77 can be a good dark matter candidate, 
and the usual seesaw mechanism does not operate because the neutrino Yukawa interactions 
can not generate masses. However, the light Majorana neutrino mass matrix can be gen- 
erated radiatively through one-loop with the help of the Yukawa interaction LiNRfi in the 
Lagrangian, which will be discussed more in detail in Sec. HI. In our model, the A^ fiavor 
symmetry is spontaneously broken by A4 triplet scalars x- From the condition of the global 
minima of the scalar potential, we can obtain a vacuum alignment of the fields x relevant 
to achieve our goal. 

The most general renormalizable scalar potential of $, 77 and x invariant under SU(2)l x 
f/(l)y X A4^ X Z2 is given as 

V = V{r]) + K($) + V{x) + V{r]^) + V{vx) + V{^x) (9) 



where 



vix) 
v{vx) 



y($x) 



+ KiV^V)3SV^V)3a + {'^5(^^^)3.(^"^^)3. + h.C.} , 

= fJ'lixx)! + >^i{xx)i{xx)i + >^2{xx)r{xx)i" + A3(xx)3,(xx)3. 
+ A^(xx)3„(xx)3„ + >^5ixx)s.Axx)3a +^iX{xx)3s +^2X{xx)3a . 

= Af (r/tr/)i($t<l>) + Xf{r]^<!>){^^7]) + hf{r]^^){7]^^) + h.c\ 

= \T{V^VUXX)1 + Ar(r/tr/)i,(xx)i" + XT{v^v)v'{xx)v 

+ Af (r/'^r/)3^(xx)3.(xx)3. + Ar(^"^^)3.(xx)3„ + Af (r/'^r/)3^(xx)3<. 

+ iTiv^vh.x 



(10) 



Here, /x^, //$, /i^, ^f , ^2'! ^1''^ and ^g'^ have a mass dimension, whereas A^ 5, A*, A^ 5, A^ 3, 
A^^ g and A*^ are all dimensionless. In ^(77$), the usual mixing term $''"77 and ^'^rjx are 
forbidden by the ^4 x Z2 symmetry. The vacuum configuration is obtained by vanishing of 
the derivative of V with respect to each component of the scalar fields $ and Xi but with 
{rji) = (z = 1, 2, 3) as follows; 



dV 



dxi 
dV 



2v^ 



dx2 
dV 



dx3 



<Xi>=-"x 



<Xt>=-"x 



<Xi>=-"x 



2 



(2A^ - A| + 4Af)(t;J^ + v^) + 2(A^ + A|)^J J 
(2A? - A| + 4Af)« + t;^3) + 2(A^ + A|)t;jJ 



^xxj^A^'^ + ZiJ 
+ 6Cft;x2t^x3 = , 

+ 6Cf^xi^X3 = 

= 2^;,3{4A*'^ + /ij + (2Af - A| + 4Af )(^;J^ + v^) + 2(A^ + A^). 

+ 6^1X1 ^x. = . (11) 



'Ms 



From those equations, we can get^ 



(Xi) = v^ 



-^x 



:A*x 



2(A? + A^) 



7^ , (X2) = (X3) = 



(12) 



^ There exists another nontrivial solution (x) 
But, it is not desirable for our purpose. 



«x(l'l'l) ^ith Wx 



-3Ji^±^9gf ~8(M^+t-|A'^x)(3A?+4A^) 



4(3Af+4AJ) 



where v^ is real. Requiring vanishing of the derivative of V with respect to $, 



dV 



d^^ 



<ip">=v^ 



2t;${24A* + 4 + X^Hvl, + vl, + vl,)} = 



(13) 



and inserting the results given by Eq. (TT2l) . we obtain electroweak VEV, 



u = f $ 



-fil - v^X'^x 



2A* 



(14) 



In our scenario, we assume that v^ is larger than w$. 

After the breaking of the flavor and electroweak symmetries, the vacuum alignment in 
Eq. fll2p leads to the right-handed Majorana neutrino mass matrix expressed as 



Mr = M 



( 



\ 



1 











1 


Ke'^ 





Ke'^ 


1 ) 



(15) 



where k = \X^v^/M\. In addition, the charged lepton sector has a diagonal mass matrix 
nii = V Dia.g.{ye,y^,yr)- We note that the vacuum alignment in Eq. ( [T2l) implies that the 
A4 symmetry is spontaneously broken to its residual symmetry Z2 in the heavy neutrino 
sector since (1,0,0) is invariant under the generator S in Eq. (E]). 

After the scalar fields get VEVs, the Yukawa interactions in Eq. ([7]) and the charged 
gauge interactions in a weak eigenstate basis can be written as 

- £ = h^MnNn + hm^R + WY.vNn + -^W-h'y'^UL + h.c , (16) 

where fj = Dia.g.{fji,fj2, fj^)- One can easily see that the neutrino Yukawa matrix is given as 
follows; 



^ y^ ^ 



Y, = V3 



( 



Ul 



with LL 



1 

7! 



V 



1 1 i\ 

\ bj bJ^ 



(17) 



yl 

For our convenience, let us take the basis where heavy Majorana neutrino and charged lepton 
mass matrices are diagonal. Rotating the basis 



^R -> UiNn , 



(18) 



the right-handed Majorana mass matrix Mr becomes real and diagonal by a unitary matrix 
Ur, 

Mr = UIMrUr = MDiag.(a, 1,6), (19) 



where a = a/1 + k"^ + 2/t;cos.^ and b = a/1 + k^ — 2/t;cos.^ with real and positive mass 
eigenvalues. Mi = Ma, M2 = M and M3 = M 6. The unitary matrix Ur diagonalizing Mr 
given in Eq. (TT5|l is 

1 



f/ 



R 



V2 



1 
1 



-1 



/ 



V 



1 

e'^ 



(20) 



with the phases 



tan ^ 



-Ksin^ 



and 'i/'2 = tan '^ 



J 



Ksin^ 



(21) 



1 + K COS ^ / V 1 — K COS C, ' 

The phases V'1,2 go to or vr as the magnitude of k defined in Eq. flTSj) decreases. Due to 
the rotation (fTSl) . the neutrino Yukawa matrix y^ gets modified to 

Yv = YMr , 

= P,tDiag.(K|,||/2l,||/3l)t/it/«. (22) 

Absorbing P^ into the neutrino field ul and then transforming ii — )■ P*iL , ^r — ;■ Pu^r, 
we can make P,^ disappeared in Y"^ as well as the Lagrangian Eq.(I7]). Then, the neutrino 
fields ul in the weak basis are simply transformed into the mass basis by the lepton mixing 
matrix, f/pMNS? so-called PMNS mixing matrix. 

The lepton mixing matrix t/pMNS can be written in terms of three mixing angles and three 
CP-odd phases (one for the Dirac neutrino and two for the Majorana neutrino) as follows 

Q 



/ 



u, 



PMNS 



Sige-^'^C'P 



\ 



-C23S12 - S23Cl2Sl3e*'''^^ C23C12 - S23Sl2Si3e*'^c;-P 



V 



,«<5cp 



Q. 



(23) 



/ 



C13C12 C13S12 

S23C13 

S23S12 - C23Ci2Si3e""'-^ -S23C12 - C23Sl2Si3e""'-^ C23C13 

where Sij = sin6'jj and Cij = cos%, and Qu = Diag.(e~*'^^/^, e"*'^^/^, 1). Here, we notice that 
the origin of the CP phases in L^pmns is the CP phases ipi,ip2 (or ^) originally coming from 
Mr as can be seen by comparing Eqs. (I15tl22p . Thus, we expect that there can be some 
correlation between low energy CP violation measurable from neutrino oscillations and high 
energy CP violation responsible for leptogenesis in the neutrino sector. 
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FIG. 1: One-loop generation of light neutrino masses. 



III. NEUTRINO MASSES AND MIXING ANGLES 



We now proceed to investigate the low energy neutrino observables. Due to the auxihary 
Z2 symmetry, the usual seesaw mechanism does not operate any more, and thus light neu- 
trino masses can not be generated at tree level. However, similar to the scenario presented 
in |8| , the light neutrino mass matrix can be generated through one loop diagram drawn in 
Fig. [1] thanks to the quartic scalar interactions. After electroweak symmetry breaking, the 
light neutrino masses in the flavor basis where the charged lepton mass matrix is real and 
diagonal are written as 



[m 



ujaji 



y-A<gOU^^ 



where 



f{z^ 



1 + 



167r^ 



Zj In Zj 

1 - Zi 



M, 



Mf 



m 



Vi 



(24) 



Ami 

'h 



\m 



R^ 



m 



Av^Xt' 



(25) 



with Zi = M^ /mi,. The explicit expressions for irti, are presented in the Appendix. Here, 

ffii ± A??7,^ /2 where 



iI^) 



m,ji-(mi.) is the mass of the fleld component rjj^Xrjj) and mjj 

the subscripts R and / indicate real and imaginary component, respectively. With Mji = 

Diag(Mri, Mj.2, M^s) and M^ = Mif~^{zi), the above formula Eq. flMl) can be expressed as 



2\<S>V 



v'X 



m„ 



( 



mo 



-YyM^ 1; = f/pMNS Diag.(mi,m2,m3)f/p]y[Ns 

^yl Byiy2 Byi ^ 
Byiy2 Dyl Gy2 
Byi Gy2 D j 



(26) 



9 



where nii^i = 1, 2, 3) are the hght neutrino mass eigenvalues, 1/1(2) = yim/Vs, and 

a a 

G ./(.,) + £!^ ^ !qM . (27) 

It is worthwhile to notice that in the limit of 2/2 —^ 1 the above mass matrix in Eq. fl26l) goes 
to fi — T symmetry leading to ^13 = and ^23 = — 7r/4. Moreover, in the limit oi yi,y2 ^ 1 
the above mass matrix gives TBM angles and mass eigenvalues, respectively, 

^13 = 0, ^23 = ~ , Ou = sin^i (^ 

mi = Smo^^e'^' , 7712 = 3mo/(^2) , ^3 = 3mo^^e'^^'+^^ , (28) 

a 

indicating that mass eigenvalues are divorced from mixing angles. However, recent neutrino 
data including the observations of non-zero ^13 requires deviations of yi^2 from unit. 

Now, let us show how deviations of 2/1^2 from unit are responsible for non- vanishing 6*13, 
and they are related with neutrino mass eigenvalues. To separately obtain real values for 
the neutrino mixing angles and masses, we diagonalize the hermitian matrix mi^ml with rriij 
given by Eq. ( I2B]) . 

/ Ayf ym (^ - ^'-^) y. (^ - ^^^^ \ 

^PMNS Diag.(mi,m2,m^) f/i^MNS ' (29) 



mi^ml = ml 



10 



where A, D, F, G, K, P, Q, R and 5* are real : 



A = (l + Ayf + yl 



n^i) 



+ (1 + yi + yi)f\z^) - 2(1 - 2yi + yi 



/(^l)/(^2) 



COSlfj 



1 ) 



[i + ^yl + ylY 



— + 4(1 + yl + yi)f\z,)+ A{l-2yl + yl)li^^^^^cosi;,, 



2 I „.2^/ l^iy I ^/i I „.2 , „.2^^2^^ \ I An <-)„.2 
i^ = 6(1 - 2/2 



/(^a) //(^i) 



COS'?/'l2 + 2/(2:2) COS ?/;2 



G = 9(1+2/2 



62 



D 
P 



[l-yl)^(^-^simPu-2f{z,)siml;,] , 



6 V a 

,2 I „,2\ 



:i+%,^+2/2^)^ 



a^ 



Q = 3(1-2/2)^7— 



a 



+ 2(1 + yl + yDfiz,) - (1 - 2yl + 2/^)^^^^^^^^ 



COS?/'i2 - f{z2)cOS1p2 



COS-?/;! 



D n O 2 , 2n/(^i)/(^2) . , 

P = (1-22/1 +2/2) ^ sm^i 



a 



^ 






sin?/^i2 + /(2;2)sinV^s 



(30) 



b \ a 

with -j/jij = t/jj — ^j. To see how neutrino mass matrix given by Eq. (l26p can lead to the 
deviations of neutrino mixing angles from their TBM values, we first introduce three small 
quantities e^, (2 = 1-3) which are responsible for the deviations of the Oj^ from their TBM 
values ; 



9' 



23 



TT 



6*13 = £2 , 6^12 = sin M — ] + €3 . 



(31) 



Then, the PMNS mixing matrix keeping unitarity up to order of e^ can be written as 



f/i 



PMNS 



/ 



v 



V2- 



es 



1+€3V^ 



£26 



-iSr 



\ 



V3 Vs 

\/6 v^ Ve Ve V2 

-l+(;i+e3\/2 _ _e2_pi5cp \/2-t3 — \/2ti _ _e2_„iScp l+ei 



V6 



V3' 



Ve 



V6^ 



g, + 0{e^) . (32) 



V2 / 



The small deviation ei from maximality of atmospheric mixing angle is expressed in terms 
of the parameters in Eq. fl5U|) as 

R{l + y2)-S{l-y2) 



tane 



P(l- 2/2) -5(1 +2/2) ■ 
The reactor angle 613 and Dirac-CP phase 6cp are expressed as 

yiM 



(33) 



tan 2^13 
tan 5cp = 3 



v^(e - A) ' 

(P-g)2 + 2/|(P + g)2 
{P + Q){R-S)-yi{P-Q){R+S) 



(34) 



11 



where 

n = {1- y2)P + (1 + y2)Q + ei{(l + y2)P + (1 - y2)Q 

- ?>i[r{1 - y2) - S{1 + y2) + ei (i?(l + 1/2) - 5(1 - 1/2)) } , 

= \[{F + G-k)(^-^ + e,{l-yl)^-y,{F-G)^ . (35) 

In the hinit of yi,y2 ~^ Ij the parameters Q,R,S,ei go to zero, which in turn leads to 
6*13 — )■ and 6cp — )■ as expected. Finally, the solar mixing angle is given as 

tan 2^12 - r- ^^ , (36) 

where the parameters \E'i, \l/2 and Z with |ej| <^ 1 are given as 

*i ^ ^-^' ^^P(i + y2) + g(i-y2)-ei{p(i-y2) + g(i + y2)} , 

^2 ^ ^ (1 + yl) + -^y2 - ei ^ (1 - yl) ■ (37) 

Note that in Eq. (IMD the condition P(l + 1/2) + <5(1 - ^2) > |ei {P(l - ^2) + Q{1 + 2/2)} | 
should be satisfied, in order for the solar mixing angle 612 to be lie in the allowed region 
from the experimental data given in Eq. (j3]). The squared-mass eigenvalues of three light 
neutrinos are given by 

ml c^ mil s?2^i + c?2^2 + Y/K sill 2^12 
2 2 I ^ , ^2!^! 



mi^ml^Q + ^j . (38) 

We see from Eq. (1371) that the deviation 63 from tri-maximality of solar mixing angle is 
roughly expressed as 

yi3V2Zml r- 

""^^"^a;;^ — ^ ^ ^ 

In the limit of |ei| ^ 1, the solar and atmospheric mass-squared differences are roughly 

given in a good approximation by 

2 

777 r ~ ~ ~ ~ ~ ~ 

Am^i^m^-m? ^ _||(F + G - A')(l + y^') + 2y2(P - G) - SA 

+ 162/i(p(l + y2) + g(l-y2))} , 
AmL^ ^ m^ - m? ^ "^{ ^ ^ \- ^\ \ + yl) - y2{F - G) - 2A 

- y,(^P(l+y,)+Q(l-y,)^^ . (40) 

12 
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FIG. 2: Plots for Case (i) displaying the reactor mixing angle ^13 versus the ratio yi/y2 (upper left 
panel), and the Jarlskog invariant Jcp versus the reactor angle ^13 (upper right panel). Allowed 
values for the atmospheric mixing angle 623 (lower left panel) and the solar mixing angle 9i2 (lower 
right panel) versus the mixing angle ^13, respectively. The thick line corresponds to ^13 = 8.68° 
which is the best-fit value of Eq. ([2]) including the Daya Bay result. And the horizontal and vertical 
dotted lines in both plots indicate the upper and lower bounds on ^13 given in Eq. ([2]) at Scr 



Here we note that the parameter Mri in Eq. (126|) can be simplified in the following limiting 
cases as 



M„~ 



Mi [In Zi-1] \ for Zi > 1 



fhlM-\ 



for Zi ^ 1 
for Zi <^ 1 . 



(41) 
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IV. NUMERICAL RESULTS 

As is well known, the observed hierarchy |Am^^.jj^| ^ ^'^soi > ^ leads to two possible 
neutrino mass spectrum: (i) rrii < m2 < m^ (normal mass spectrum), and (ii) m^ < rrii < 
7712 (inverted mass spectrum). Since there are many unknown parameters such as masses 
of heavy Majorana neutrinos and scalar fields r//?,?]/, we consider a particular parameter 
set for those parameters and show how the measured values of the mixing angle 6*13 can be 
accommodated in our model while keeping the other neutrino parameters such as solar and 
atmospheric mixing angles and mass-squared differences are satisfied with the current data. 

The mass matrix in Eq. 0261) contains 10 free parameters : X^^,M,y^, Zi,Z2,z^ and 
2/1)2/2; 'C) f^- The combination of the first three of them, {A3 '', M, 1/3 }, leads to the overall neu- 
trino scale parameter mo. As shown above, the elements of the mass matrix in Eq. fl26l) are 
expressed in terms of measurable neutrino parameters, 9^, 9rs, 623, rrii 2,3, Sep, ^1,2- Among 
them, three mixing angles and two mass squared differences are measured. For numerical 



analysis [ll| , we need to fix some parameters by hand since there are too many model param- 
eters to be predicted. As an example, we take a case Mf = fn^^, M2 = l.Srh'^^ M^ = l.Sm^^, 
and fix the overall seesaw scale M to be 1 TeV. Then, the parameters tuq, yi,y2,i^,^ can be 
determined from the experimental results of three mixing angles and two mass squared dif- 
ferences. In addition, the CP phases 6cp, V5i,2 can be predicted after determining the model 
parameters. Depending on the values of the model parameters, there exist two possibilities 
for the light neutrino spectrum, one is normal mass hierarchy and the other is inverted 
hierarchy. In the following, we discuss the two cases separately. 

(i) normal hierarchy of light neutrino 

Based on the formulae for the neutrino mixing angles and masses, we numerically scan 
the parameters mo,yi,y2, n,^ and then pick up the values of those five parameters which 
are consistent with the experimental data given at 3a in Eq. ([3]). For the mixing angle ^13, 
we a bit widely allow its value from 5° to 15° instead of its experimental values at 3a^ . In 



^ Note that very small mixing angle 6*13 less then 1° can be achieved in the case that ?/i — > 1 or sin^'i — > 
converges more faster than 2/2 ^^ 1- 
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FIG. 3: Predictions for the Dirac CP phase 5cp versus ^13 (left panel) and the effective mass of 
neutrinoless double beta decay |?7iee| versus the mixing angle ^13 (right panel) for Case (i) . The 
thick and dotted lines correspond to ^13 = 8.68° which is the best-fit value and the 3cr bounds 
given in Eq. ([2]) including the Daya Bay result, respectively. 

such a way, we can obtain the allowed regions of the parameters given by 



1.40 <K< 2.38 
190° < ^ < 211° 



0.44 < yi < 0.89 



0.60 < 2/2 < 0.84 and 1.1 < y2 < 1-89 



0.23 < 



10 



^ < 0.46 . 



-9 



(42) 



We found that normal mass ordering of light neutrino can be achieved when Mi < M2 < M3 
or M2 < Ml < M3 are satisfied for the parameter spaces given above. In the left up- 
per panel of Fig. [21 the data points indicate how the mixing angle ^13 is determined in 
terms of the ratio yxjyi- The result shows that the upper limit of 2/1/2/2 is 0.86, and the 
measured value of 6'i3 from the Daya Bay and RENO can be achieved for two regions, 
0.40 < 2/1/2/2 < 0.57 and 0.67 < 2/1/2/2 < 0.82. To see how the parameters are correlated 
with low energy CP violation measurable through neutrino oscillations, we consider the lep- 
tonic CP violation parameter defined by the Jarlskog invariant Jcp = Im[?7eif/^2f^e2^^i] = 
I sin 26*12 sin 2^23 sin 26*13 cos 6*13 sin 6cp |12| which can be described in terms of the elements 



h = myml 13| : 



. _ Im{/li2/l23/^3l} 

Am2iAm3j^Am|2 



(43) 



The behavior of Jcp is plotted in the right upper panel of Fig. |2]as a function of 6*13. We 
see that the value of \Jcp\ lies between and 0.034 for the measured value of ^13. In 
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our model, since Im{/ii2/i23^3i} is proportional to 1 — y|, the leptonic CP violation Jcp 
goes to zero in the limit of 1/2 — ^ 1- However, 2/2 = 1 is not allowed in our analysis, and 
thus Jcp = indicates that there exists some cancelation among the terms composed of 
sin?/^i2,sin(?/;i + ^2), sin(2?/'i - ^2) and sin?/'2 multiplies by 2/1,2, f{zi)/a, f{z2), and f\z^)/h 
even if CP phases il)i^2 are non zero. In the lower panel of Fig. [21 the data points indicate how 
the values of ^13 depend on 6^12 and ^23 in the allowed regions given by Eq. (|3]). We see that 
the measured values of 61^ can be achieved for two separate regions of 6^23 : 38.6° < ^23 ^ 43° 
and 47° < ^23 ^ 53.1°, which indicates that the parameter set strongly prefers deviations 
from maximal mixing for the atmospheric neutrino oscillation. From the right lower panel of 
Fig. 121 we see that predictions of ^13 does not strongly depend on 612 for the allowed region. 
We see from the figures that 9i^ for the normal hierarchy prefers rather large values more 
than 5 degrees. 

We also see from Fig. |2]that small deviations for 6'23 prefer to large value of ^13 in normal 
hierarchical case. This can be understood by considering two relations given in Eq. (21) and 
Eq. (26). The phases ipi^2 go to or tt as the magnitude of k defined in Eq.(15) decreases, and 
in the case of ?/2 = 1 the neutrino mass matrix indicates directly 6'i3 = and ^23 = — vr/4. 
However, deviation of 2/2 from one can be associated with deviation from maximality of 
atmospheric mixing angle by the following relation, 

l + y^\ {l-2yl + yl)^mi;,l^^^^-{l-y2f^(^^mij,2 + f{z2)^ini;2) 
tanei — ' 1 



1 - ^2; (1 _ 2yl + yl) sin^.lMli^ _ (1 + y^)2li^ (^/(£i) sj^^^^ + /(^2) sin ^2 

This formular for the parameter ei is relevant only when ^2 7^ 1- In the case of ^2 — ^ 1 
while yi ^ 1 and saiipi 7^ 0, we see from the above equation that the value of ^23 (or ^i) 
can be large but restricted by experimental data. Then, due to Eq .(26) and Vt in Eq. (35), 
the value of ^13 gets smaller as ^2 — ^ 1- On the other hand, when 2/2 is much deviated from 
1, two cases for 6'23 (or ei) are possible. One is that rather smaller values of ^23 (or ei) are 
preferred as the value of k (or sin-?/;! -> and sin'?/'2 -> 0) decreases, and the other is that 
the combination of two parts in numerator of the above equation can lead to wide ranges 
of 6*23 (or ei). However, when 2/2 ~ 1 which makes the above equation irrelevant, the value 
of 6*13 goes to 0° (numerically < 1°), and the value of ^23 can approach 45° (or ei — )■ 0) for 
1/1 — )■ 1 or sin'?/'i -^ converge more faster than ?/2 — )■ 1. We have neglected this case in our 
paper. 
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By using the conventional parametrization of the PMNS matrix 10| and Eq. (132|) . one 
can deduce a expression for Dirac CP phase 5cp given by 

dcp = - arg '- . (44) 

\ S12S23 / 

Moreover, we can straightforwardly obtain the effective neutrino mass \mee\ which is asso- 
ciated with the amplitude for neutrinoless double beta decay : 



irrip 



5^(t/pMNs)ei"^i| , (45) 



where t/pMNS is given in Eq. (15^ . The left panel of Fig. [3] shows that 6cp is predicted to 
be 0° < 5c p < 60°, 120° < 5c p < 240° and 300° < 5c p < 360° for the measured values of 
^13 at 3a. In the right panel of Fig. [31 we plot the prediction of the effective neutrino mass 
|mee| as a function of ^13, which lies between 0.014 and 0.021 in the region of the measured 
values of 6^13 at 3a. 

(ii) inverted hierarchy of light neutrino 

Now let us turn to the inverted hierarchical case. Similar to case (i), scanning the 
parameters mo, yi,y2, k,,^ based on the formulae for the neutrino mixing angles and masses 
and taking the experimental data given at 3a in Eq. ([3]) as constraints, we can obtain the 
allowed regions of model parameters given by 

1.30 < K< 1.56 , 209° < ^ < 222° , 0.27 < ^Mr < 0.45 , 

0.79 <yi< 0.88 , f 1.12 < y, < 1.24 , 

and < (46) 

0.60 < 2/2 < 0.79 , [ 1.28 < 2/2 < 1.5 . 

We found that this case is achieved when Mi < M2 < M3 is satisfied. For those parameter 
regions, we in turn investigate how the mixing angle ^13 depends on other parameters and 
whether CP violation is realized. In the left upper panel of Fig. HI the data points indicate 
how the mixing angle ^13 is determined in terms of the ratio yi/y2- We see that the measured 
value of ^13 in 3a including the Daya Bay experiment in Eq. ([2]), can be achieved for two 
separate regions, 0.82 < yi/y2 < 0.88 and 1.12 < yi/y2 ^ 1.3. We plot Jcp vs. 613 
in the right upper panel of Fig. H For 5.9° < ^13 < 9.5°, \Jcp\ ^ 0.018 ~ 0.036 and 
—0.02 ~ —0.034, which indicates CP violation in the leptonic sector. 
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FIG. 4: Same as Fig.[2]for Case (ii). 

In the lower panel of Fig. IH the data points show how 6*13 is determined in the allowed 
regions of 612 and ^23 given by Eq. (|3]). We see that the narrowed regions of the atmospheric 
mixing angle ^23, 38.6° < 6^23 < 40.5° and 49.5° < ^23 ^ 53.1° are preferred, which indicates 
that the parameter set disfavors maximal mixing for the atmospheric mixing angles. From 
the lower right panel of Fig. HI we see that determination of 9i^ does not strongly depend on 
9i2 for the allowed region. We see from the figures that contrary to the case (i), ^13 for the 
inverted hierarchy prefers rather lower values less than 9.5 degrees. The left panel of Fig. O 
shows that 5cp is predicted to be around 70°, 100°, 160°, 250° and 290°. In the right panel 
of Fig. [5l the value of \mee\ is predicted as a function of 6*13 and we see that |mee|[eV] lies 
between 0.038 and 0.049 in the allowed region of ^13. 
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FIG. 5: Same as Fig. [3] for Case (ii) 



V. CONCLUSION 



Motivated by recent observations of non-zero 613 from the Daya Bay and RENO experi- 
ments, in this paper, we have proposed a neutrino model with ^4 symmetry and shown how 
deviations from the TBM mixing indicated by the current neutrino data including the Daya 
Bay result can be accounted for. In addition to the leptons and the Higgs scalar of the SM, 
our model contains three right handed heavy Majorana neutrinos and several scalar fields 
which are electroweak singlets and demanded to construct desirable forms of the letponic 
mass matrices. To have a good dark matter candidate, we imposed auxiliary Z2 symmetry, 
and thus light neutrino masses at tree level are absent in our model. However, the light 
neutrino masses can be generated through loop diagram, and we have shown how the light 
neutrino mass matrix can be diagonalized by the PMNS mixing matrix whose entries are 
determined by the current neutrino data including the Daya Bay result. In our model, 
the origin of the deviations from TBM mixing is non-degenerate neutrino Yukawa coupling 
constants among three generations. Also, unremovable CP phases in the neutrino Yukawa 
matrix are the origin of the low energy CP violation measurable from neutrino oscillation 
as well as high energy CP violation. We have discussed some implication on leptonic CP 
violation. 
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Appendix A: The Higgs mass 

Our model contains four Higgs doublets and three Higgs singlets. Here, we present 
the masses of physical scalar bosons, where the standard Higgs h' is mixed with Xqj, 
not with h'^^A'^. For simplicity, we assume that CP is conserved in the scalar poten- 
tial, and then the coupling A3 is real and the term ^2^{i]'^f])3aX is neglected in the 
Higgs potential given in Eq. flTOj) . The neutral Higgs boson mass matrix in the basis of 
(^'j Xoi7 Xo27 Xo37 ^i! ^2) ^3) ^1) ^25 ^3) is block diagonalized due to Z2 symmetry and CP 
conservation, which is given by 
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(Al) 



/ 



where the primed particles are not mass eigenstates, and mass parameters are given as 



m,, 



Vfl 

m 



xi 
2 



m\ 



VTiu 



VTiv 



m , 



m, 



m 



h'^K 



AXW 



'$ : 



^h'x[ = 2t^<J>^x^*'' 



Avl{X^ + Af ) , ml 



^(3A^ 



4A^: 



X^ -^ ■ ^' ' X2(3) X^ ^ ■ o, . 

vUxf + Xf + 2Af ) + /ij + vliXr + 2Re[Af ]) , 
vUxf + Xf - 2Af ) + ^; + vliXr + 2Re[Af ]) 



^kx's 



3^xe? 



+ xf + 2Af ) + /ij + vliXr - Re[Af ] - v^lm[A^^]) , 



4(Af 

vUxf + Xf + 2Af ) + /ij + v^iXr - Re[Af ] + y3Im[Af ]) , 



liXf + xf - 2Af ) + /i; + vliXr - Re[Af ] - v^Im[Af ]) 
(Af + Xf - 2Af ) + /ij + vliXr - Re[Af ] + v^Im[Al^]) 



m 



— o, cvx 



2 c^i 



(A2) 
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Since the matrix in Eq. f lAip is block diagonalized, it is easy to obtain the mass spectrum 
given as follows; 



^l = 


1/ 2 


+ mj, -y/(m2, 


-^l,y + Aiml^,f 


<- 


+ mj, + ^{ml 


-mJ,)^ + 4K,^,P 


K2 -- 


= "^k- 


Kw, ' ^?3 = 


X2 X2X3 ' 


rnl, = 


= ^l[ . 


"^1 = ^\ , 





m 



h2 



vl{\t + 2Af ) + /i? + ^J(Ar - Re[Af ]) - v^sJ'i{v^^^{\t\Y + (e? 



J'7XN2 



A2 

2 
'/»3 



4(A?* - 2Af ) + /ij + ^J(Ar - Re[Af ]) - v^p{v^^^{\f\f + (O^ , 
4(A?* + 2Af ) + /ij + vl{\r - Re[Af ]) + v^^j'^{y^^^\\f\f^{iTY , 
m% = 4(A?2* - 2Af ) + /ij + t;J(Af - Re[Af ]) + v^^ 3{v^Re[Xf]y + (^D^ , (A3) 
where A^2 — -^i + -^2 • Note here that the unprimed particles denote mass eigenstates. 



And the charged Higgs boson mass matrix in the basis of (j?! ,^72 jW ) is given as 



/ 



m 



charged 



m\ 

m2± 

V2 



\ 



(A4) 



\ m^i / 



where 



m 

Vi 

rn?± 
V2 



± = /i^ + 4Ar + i;^(Af + 2Re[Af]) 



/i; + vlXf + vl (Af - Re[Af ] - v^Im[A 
m2± = fil + vlXf + vl (Af - Re[Af ] + V3lm[X 



»?xi 
2 J 



(A5) 



Using ml., m\^ in Eq. ( JA2I) and Eq. ( lA4p . the expressions for ?72^. appeared in Eq. (I24 



are 






^l + 4Af + t;^ (Af + 2Re[Af ]) = m^ + vlXf 



V. 



2 \V^ 



Vl 



A?2* + f^l + vliXr - Re[Af ] - v^Im[Af ]) = m\ + ^A 



V2 



m. 



V3 



V. 



2 \ri^ 



Kt + /i' + vliXT - Re[Af ] + v^Im[Af ]) 



m\ + v'ixT 



2 \V^ 
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(A6) 
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